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FREE-FLIGHT INVESTIGATION 

AT SUPERSONIC SPEEDS OF THE STABILITY AND DRAG OF A 
79° CLIPPED DELTA BOOST -GLIDE CONFIGURATION INCLUDING 
AN ANALOG STUDY OF COUPLED MOTIONS 
DURING THE FLIGHT* 

By Sherwood Hoffman and Willard S. Blanchard, Jr. 

SUMMARY 


A free-flight investigation was conducted at Mach numbers between 
1.2 and 3. A to determine the aerodynamic characteristics and motions 
of a rocket -boosted model of a simplified boost-glide hypersonic con- 
figuration. The model had a 79° clipped delta planform, sharp leading 
edges, a blunt base-, and symmetry in two planes. The mass moment of 
inertia in roll was about 4 percent of the value in pitch or yaw. 

The model was both statically and dynamically stable near Mach 
number 3 at small angles of attack and slideslip. A disturbance in 
pitch to an angle of attack of about 12° resulted in roll oscillations 
and reversals, and coupling of the longitudinal and lateral forces. 

An analog study, in which the equations of motion for five degrees of 
freedom were used, generally simulated the magnitudes and frequencies 
of the angular motions and showed that the rate of change of effective 
dihedral derivative with angle of attack was a necessary derivative for 
simulating the coupled motions at the altitude and velocity of the test. 

Linearized theory gave reasonable predictions of the pressure drag 
and of the static and dynamic stability derivatives for oscillations of 
small amplitude. 


*Title, Unclassified. 
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The National Aeronautics and Space Administration is devoting con- 
siderable effort to the design of boost-glide vehicles for flight into 
space, reentry, and landing. Previous studies, such as references 1 
and 2, indicate that the reentry phase of such a mission is feasible 
from both aerodynamic -heating and loading aspects with properly designed 
lifting configurations. Wind-tunnel tests have been made to determine 
the aerodynamic characteristics of promising boost-glide shapes, some 
examples of which may be found in references 3 to 5 for lifting bodies, 
references 6 to 8 for all-wing vehicles, and reference 9 for wing -body 
combinations. The results presented in the references show that low- 
aspect-ratio configurations may experience significant rolling moments 
due to aerodynamic coupling of their longitudinal and lateral forces 
at angles of attack or sideslip. The resulting rolling motions may 
become pronounced because the mass is concentrated along the longitudi- 
nal axis and therefore the inertia in roll is very small compared with 
the inertias in pitch and yaw. Information on the motions and on the 
cross-coupling derivatives involved would be valuable to the designers 
of boost-glide configurations, especially with regard to controls. 

A rocket-boosted model of a simplified hypersonic glider configura- 
tion has been flight-tested to determine its aerodynamic characteristics 
and motions during free flight at supersonic speeds. The model was 
representative of an all-wing boost-glide vehicle in that its leading 
edge was highly swept and the maximum thickness was at its blunt base. 

The configuration differed from some proposed boost-glide configurations, 
since the leading edge was fairly sharp instead of rounded and because 
it had symmetry in two planes. 

The model was flight-tested at the NASA Wallops Station, through a 
range of Mach numbers decelerating from J.k to about 1.2 with corre- 
sponding Reynolds numbers, based on mean aerodynamic chord, ranging 
from 30 x 106 to about 2.5 x 10°. Pulse rockets disturbed the model 
in the pitch plane during the flight and a ten-channel telemeter trans- 
mitted continuous acceleration and pressure data to ground receiving 
stations. The analysis includes the determination of stability deriva- 
tives from small-amplitude oscillations, comparisons with predictions 
from linearized theory, and an analog simulation study of the coupled 
motions during large -amplitude oscillations. A small model also was 
tested from a helium gun to determine the transonic drag of the 
configuration. 
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SYMBOLS 


The "basic data are presented with respect to the body-axis system 
shown in figure 1. The axes originated at the center-of -gravity posi- 
tion located at the 55 .1-percent longitudinal station on the longitudi- 
nal axis of the rocket model, and corresponded to the J>1. 8-percent 
station of the mean aerodynamic chord. 


a L,cg 

a L,Hi 


a- 


L,Lo 

^cg 

a N,n 

a N,t 


%cg 

b 


longitudinal acceleration of center of gravity, g units 

longitudinal accelerometer reading, high range, positive 
in positive X direction, g units 

longitudinal accelerometer reading, low range, positive 
in positive X direction, g units 

normal acceleration of center of gravity, g units 

normal accelerometer reading in nose section, positive 
in negative Z direction, g units 

normal accelerometer reading in tail section, positive in 
negative Z direction, g units 

lateral accelerometer reading, positive in positive Y 
direction, g units 

lateral acceleration of center of gravity, g units 

span, ft 

mean aerodynamic chord, ft 


°D 

C D,b 

C D,f 

C D,p 


drag coefficient, 


Drag 
q S 


base-drag coefficient, based on S, ( 5s — ] x — x 144 

\ %, ) S 

friction-drag coefficient, based on S 
pressure-drag coefficient, based on S 


CONFIDENTIAL 




* • • • • • • • 

• • • • 
• • • • • • 

> • • • • • 
• •• ••• • • 


• • • 
• • • 


CONFIDENTIAL 


lift coefficient, 


rolling -moment coefficient, 


Rolling moment 


pitching -moment coefficient, Pitching mo m en t 

<3co S5 


normal-force coefficient , 


Normal force 


yaving -moment coefficient, 


Yawing moment 


axial-force coefficient, positive in positive X direction, 

Axial force 
qooS 


lateral-force coefficient, 


Lateral force 


lift-force coefficient derivative, , P er radian 


effective dihedral derivative, — -, per radian 

dp 


rate of change of effective dihedral derivative with angle 


of attack, it P er square radian 

6a 


damping -in -roll derivative, — — , per radian 

6 ^ 


rate of change of rolling-moment coefficient with yawing 

6 Ci 

angular velocity factor, — — , per radian 
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m. 


a 


dC. 


static stability derivative, — per radian 

da 


Sx 


rate of change of pitching-moment coefficient with rate 

dCm 

of change of angle-of -attack factor, — — , per radian 

a ac 
d 27 


mn 


rate of change of pitching -moment coefficient with angle 
dC m 

of sideslip, — per radian 

^3 


m r 


rate of change of pitching -moment coefficient with 

dCm 

pitching angular velocity factor, — — , per radian 

2V 


"N 


a 


c3Cm 

normal-force coefficient derivative, — per radian 

da 


' n P 


dc n 

directional stability derivative, — - per radian 

dp 


C n . 

n p 


rate of change of yawing -moment coefficient with rate of 

dC n 

change of angle-of-sideslip factor, — , per radian 


n-r 


rate of change of yawing-moment coefficient with rolling 

dC 

angular velocity factor, — h_., per radian 

d 2V 


v n„ 


rate of change of yawing -moment coefficient with yawing 

dC n 

angular velocity factor, — — , per radian 

d ^ 

2V 
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lateral-force coefficient derivative, 


dC v 

— i, per radian 
dp 


rate of change of side-force coefficient with angular 

dCv 

velocity factor in yaw, — • — , per radian 

a ik 

2V 

acceleration due to gravity, 52.2 ft/sec2 
altitude, ft 

moments of inertia in roll, pitch, and yaw about the X, 
Y, and Z axes, respectively, slug-ft 2 

free- stream Mach number 

mass, W/g, slugs 

period of short-period oscillation, sec 

angular velocity in roll, pitch, and yaw, radians/sec 

angular accelerations in roll, pitch, and yaw, 
radians/sec 2 

base pressure, lb/sq in. 

free-stream static pressure, lb/sq in. 
free-stream dynamic pressure, lb/sq ft 
Reynolds number, based on c 
planform area, sq ft 
base area (including fins), sq ft 

time for a transient oscillation to damp to one -half 
amplitude, sec 

time, sec 

free-stream velocity, ft/sec 
weight, lb 
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X,Y,Z 

x,y,z 


a 

a 

3 

3 

7 


body-axis system, ft 

distances from center of gravity, measured along X, Y, 
and Z axes, ft 

angle of attack, radians 

da, 

rate of change of angle of attack with time, - — , 

dt 

radians/sec 

angle of sideslip, radians, or Y M^ - 1 

A r\ 

rate of change of angle of sideslip with time, — , 

dt 

radians/sec 

angle between flight path and horizontal, deg 
MODELS AND INSTRUMENTATION 


Drawings of the two models tested are presented in figure 2. The 
larger model (fig. 2(a)) was designed to be boosted by rockets whereas 
the smaller one (fig. 2(b)) was scaled down to permit it to be propelled 
from a 6-inch helium gun. Longitudinal distances shown are measured 
from the "theoretical" nose tip, which is 0.66 inch ahead of the actual 
nose. Photographs of the models are shown in figures 5 and 4. Physical 
characteristics of the rocket model are presented in table I. The dimen- 
sions of the small helium-gun model were 0.181 of those of the rocket 
model. The weight of the small model was 1.128 pounds and the center of 
gravity was located at the 42.5-percent station measured from the theo- 
retical tip. The thin, pointed nose tip was strengthened with a small 
nose sting as is shown in figures 2(b) and 4. 

The configuration had mirror symmetry in the horizontal and verti- 
cal planes, sharp leading edges, and a blunt trailing edge. The planform 
was a 78*87° clipped delta wing with streamwise tips, aspect ratio of 
0.54, and taper ratio of 0.191. The side view consisted of a 4° half- 
angle triangular forebody, a flat afterbody, and a maximum thickness 
of 7*65 percent of the total length. All the surfaces were planar. The 
resulting body cross sections normal to the longitudinal axis were 
diamond-shaped from the nose to about 55 percent of the body length, 
hexagonal from the 55-pencent station to about the 80-percent station, 
and octagonal from the 80-percent station to the trailing edge or base. 
The fins, which were located in the vertical plane of symmetry, had a 
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sweepback angle of 76° along the leading edge, a taper ratio of 0 . 33 * 
an unswept trailing edge, and 7-67° total-angle wedge airfoil sections. 

A ten-channel telemeter was installed in the rocket model and 
included the following instruments: high-range and low -range longi- 

tudinal accelerometers, normal accelerometers located near the center 
of gravity and in the rear of the model, a transverse accelerometer, 
angular accelerometers measuring pitching and yawing accelerations, roll- 
rate gyro, total pressure pickup and base pressure pickup. The base pres- 
sure was averaged over the semispan of the blunt trailing edge with a 
manifolded tube as is shown in. figures 2 (a) and 3 (b)- The locations of 
the instruments with respect to the center of gravity are given in 
table II. 


TESTS 


A photograph showing the rocket model and booster on the launcher 
at the NASA Wallops Station is presented as figure 5. The booster con- 
sisted of an Honest John for the first stage and a Nike for the second 
stage. A small rocket motor was installed in the base of the model (see 
fig. 2(a)) and was programed to separate the model from the second stage 
after Nike burnout. The calculated relative deceleration of the burned- 
out Nike and separated model, based on drag-to-weight ratios, was very 
small and indicated the possibility of collision after separation. Two 

3i-inch aircraft rocket motors were therefore employed to increase the 

k 

Nike deceleration after the model and booster were separated. These 
motors were mounted on the forward end of the Nike (fig. 5) with their 
nozzles facing upstream at a 20° cant angle with respect to the booster 
center line. The nozzles were plugged to protect the grain and igniters 
from aerodynamic heating during the booster phases of flight. The model 
was disturbed only in pitch by 27-pound-second pulse rockets (burning 
time 0.05 second) approximately 5 seconds and 15 seconds after separation 
from the booster. Calculations indicated that each pulse rocket would 
pitch the model to an angle of attack of about 12°. The locations of the 
pulse rockets are shown in figure 2(a) . 

Flight-path data were obtained by tracking the model with the 
AN/FPS -16 tracking radar, the NASA modified SCR -584 tracking radar, and 
the Reeves modified SCR -584 radar. Atmospheric conditions were obtained 
from a rawinsonde balloon that was released just prior to the test. Mach 
number and dynamic pressure were determined from the total -pressure meas- 
urements on the model and the ambient pressures and temperatures at corre- 
sponding altitudes along the flight path. 
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The range of flight for which the data were reduced varied between 
Mach number 3*4 at an altitude of 35 >000 feet to Mach number 1.2 at an 
altitude of about 68,000 feet as shown in figure 6. The uncertainties 
in the aerodynamic coefficients were inversely proportional to the free- 
stream dynamic pressure and varied from a small amount at M = 3-4 

and q^ = 4,000 lb/sq ft to a fairly large percentage of the coefficients 
near M = 1.2 and q^ = 100 lb/sq ft. The corresponding Reynolds num- 
bers, based on c, varied from 30 X 10^ to about 2.5 X 10^ as shown in 
figure 7 . 

Telemetered data, static pressure, dynamic pressure, and Mach number 
for three time intervals of the test are presented in figure 8 . The 
points were machine plotted from magnetic-tape records and show occasional 
scatter points due to noise. The time interval between 29-7 seconds and 
30.6 seconds, figure 8 (a), shows that the model was disturbed by the 
booster at separation (M ~ 3*4) and experienced small-amplitude sinusoidal 
oscillations in pitch q and yaw r. The roll rate p was essentially 
constant and less than 2 radians per second. When the first pulse rocket 
fired (near M = 2.8), the angle of attack increased rapidly and then 
the model motions became coupled. The resulting roll rates initially 
exceeded the range (about ±30 radians per second) of the gyro. The inter- 
vals from 40.8 seconds to 41.7 seconds (M ~ 2 . 3 ) and from 42.2 seconds 
to 43.1 seconds (M ~ 2.2) show the on-scale measurements after the first 
pulse. (See figs. 8 (b) and 8 (c).) The accelerations and motions are non- 
linear and coupled. The model oscillated in roll to the right and then 
experienced two roll reversals in the time interval from 40.8 seconds 
to 43.1 seconds. In general, the roll-rate oscillations following both 
pulses were unsteady and varied from positive oscillations to roll 
reversals and to negative oscillations. 

The helium-gun model covered a Mach number range from 1.2 to 0.7 
with corresponding Reynolds numbers (based on c) varying from 6 X 10° 

to about 3.5 X lO^ as shown in figure 7 . Velocity, trajectory, and 
atmospheric data were measured by the CW Doppler velocimeter , the NASA 
modified SCR-584 tracking radar, and a raw insonde balloon, respectively. 


ANALYSIS 


Total Force and Moment Data 

The total force coefficients were determined from the instantaneous 
values of translatory acceleration of the center of gravity obtained by 
the method described in appendix A. The following relationships were 
used to compute the force coefficients: 
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-Cy = a r 

x q^S L > C S 

(1) 

w 

C N = 

(2) 

p w 

Y = 7s a Y,cg 
00 

(3) 

The helium-gun model was tested at zero lift and had no instrumen- 
tation*. The drag coefficients were determined from the flight-path 
angle, free-stream conditions, and decelerations obtained by a differ- 
entiation of the velocity-time curve from the CW Doppler velocimeter 
by use of the following relationship: 

C D — (— + g sin 7 'j 

gq^S V dt J 

(4) 


The total pitching- moment , yaw ing- moment , and rolling-moment coef- 
ficients for the rocket model were computed from the instantaneous values 
of angular acceleration and angular velocity with the following 
expressions : 


Iyq 

+ ^L 

- Iz) 

pr 

^ S5 


QooSc 


iz* 

+ (il 

- ix) 

pq 



QooSt 


*x* 

+ ( J Z 

- M 

qr 






(5) 

( 6 ) 

(7) 


It should be noted that the values of the product-of-inertia terms 
were zero because the principal axis coincided with the longitudinal 
axis of the model. The values of the gyroscopic reaction terms became 
large when the model roll rates were large. Values for q and r were 
determined by integrating q and r over small time intervals and 
for p by differentiating the time-history records of p. The curves 
of q and r were assumed to be centered about zero to determine their 
magnitudes and signs through several cycles. 
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Small-Disturbance Analysis 

The small-disturbance data near M = 3.4, shortly after the model 
separated from the booster, have been analyzed separately in the pitch 
plane and in the yaw plane by assuming two degrees of freedom. As was 
noted earlier, the roll rate near M = 3.4 was very small and practi- 
cally constant. The values of C m at a = 0 and C^ at B = 0 

Ia a J3 

were computed from the average period and time-to-damp to one-half ampli- 
tude of the transient oscillations in each plane by the following 
relationships : 


kn 2 i. 




kn 2 Iy 


n p 


r. + i /q.695 \ 2 
p 2 4it 2 y T i/2 J 

1_ + 1 /0.693 V 

P 2 4« 2 \ 1/2 j 


(8) 

( 9 ) 


and the corresponding values of C 


and C y were computed from 
a X 3 



(10) 



The rotary damping derivatives were calculated as follows: 




Cm • — 

a 


_2 


(- 


586V 

q«) ST i/2 


52.2 

W 



C 


n r 



2I z /1.586V + 

b 2 ^^1/2 



(11) 


(12) 


Simulation of Coupled Motions 

A simulation study of the coupled motions near Mach number 2 was 
made using the equations of motion for five degrees of freedom. The equa- 
tion associated with the drag force was omitted since the velocity of the 
model was nearly constant through the time intervals studied. The cross- 
product inertia terms have been omitted from the equations since the 
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principal axis coincided with the body axis. The equations of motion 
were written for the body-axis system and are given in appendix B. A 
highly simplified motion study using three degrees of freedom also was 
made and the equations used are included in appendix B. 


RESULTS AND DISCUSSION 


The stability derivatives and drag coefficients referred to in the 
following sections are for the conditions of zero angle of attack and 
zero angle of sideslip, except where noted. 


Longitudinal Stability 

The variations of normal-force coefficient, pitching-moment coeffi- 
cient, and axial-force coefficients with time near Mach number 3.4 are 
presented in figures 9(a), 9(b), and 9(c). The corresponding variations 
of Cjn with Cjj are given in figure 10. The oscillations are sinusoidal 
about a trim normal force of zero (due to model symmetry), damp with time, 
and have an averaged period of 0.175 second. The average slope 
in figure 10 indicates that the static margin was -0.124 and the center 
of pressure was located at the 6l. 7-percent longitudinal station measured 
from the theoretical tip of the configuration. 

The time-to-damp to one-half amplitude was obtained by 

plotting the peak amplitudes of the transient oscillations of C N 
and ^ against time on semilog paper as is shown in figure 11. The log 
decrement of the peak amplitudes was linear and T \j 2 was 0.62 second. 

The resulting values of C m and Cw were computed to be -0.104 

nx a 

and O.837, respectively. The amplitude of the pitching oscillations, 
based on these values, varied within ±0.4°. 

The static stability derivatives are compared in figure 12 with 
those from linearized theory (refs. 10 and 11) for thin, sweptback, 
tapered wings at M = 3-4 and with wind-tunnel data (ref. 9) at M = 2.01 
and 1.41. The agreement with linearized theory is good and the variation 
of Cfj and C_ with Mach number is as would be expected for low- 

% <x 

aspect -ratio wings. 

The sum of the rotary damping derivatives in pitch C m + C^. 

q a 

at M = 3*4, as determined by equation (ll) was -0.077* T n comparison, 
linearized theory (refs. 11 and 12) predicted a value of about -0.4. 

In a free-flight test of a 67.5° sweptback arrow wing (ref. 13) the 
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damping in pitch also was lower than the theoretical value at M = 1 . 8 . 
It should be noted that these magnitudes are small and that the accuracy 
of measuring may be poor because they were obtained from the 

difference of two numbers having the same order of magnitude (eq. (ll)). 


Directional Stability 

The variations of side-force coefficient and yawing-moment coeffi- 
cient with time near M = 3.4 are presented in figures 13(a) and 13 (b) 
and the corresponding variation of C n with Cy is given in figure 14 . 
The yawing oscillations were sinusoidal, damped with time, and had a 
linear variation of the log decrement of the envelope of the transient 
oscillations (fig. 15) - The rolling moment during this time interval is 
shown in figure 13(c) to be practically zero. 


The data show an average yawing period of O.27O second, ^\/2 e( ^ ua ^ 

to 1.16 seconds, static directional stability ^AC n /AC y) of -O.65, and a 

center of pressure in yaw at the 73.8-percent longitudinal station from 

the theoretical nose tip of the configuration. The resulting values 

of C n p and Cy^ were O.O98 and -O.I5I, respectively, and the amplitude 

of the yawing oscillation varied within. ± 1 . 2 °. The static directional 

stability is compared with theory at M = 3.4 and wind-tunnel data 

(ref. 9) at M = 2.01 and M = 1 . 4 l in figure l6. The same theoretical 

methods that were used in pitch were employed with the exception that the 

derivatives for the body alone and vertical fins alone were calculated 

separately and then summed. The theoretical values of C v and CL 

13 “P 

are somewhat higher than the test values and this may be due, partly, to 
the omission of interference calculations between the fins and body. 


The sum of the measured rotary damping derivatives Cn r - C n ^ from 

equation (l 2 ) was -O.652 at M = J.k. The theoretical damping due to the 
vertical fins alone was -O.69. 


Drag 

The variations of total-drag coefficient and base-drag coefficient 
with Mach number at approximately zero angle of attack and angle of side- 
slip for the rocket model, helium-gun model, and for the geometrically 
similar wind-tunnel model (ref. 9) are presented in figure 17. Since a 
and 3 were not measured in flight, the values of for the rocket 

model (at a 3 p = 0 ) were taken as those values corresponding to 

• 1 ■ — ■ 

\ + nearly zero. The drag for the helium-gun model is an 
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average drag about the zero-attitude condition since it is probable that 
the model was disturbed slightly upon ejection from the gun and oscillated 
through small angles. 

Good agreement was obtained in Cp and Cp ^ for the rocket model 
and wind-tunnel model at supersonic speeds. The agreement shown in fig- 
ure 17(a) between data for the rocket model and for the smaller helium- 
gun model near M = 1.2 is poor; and the general level of the drag curve 
from the smaller model appears to be lower in comparison. Part of this 
difference may be explained by the different test Reynolds numbers and 
their effect on the friction drag. Computed Cp ^ curves using average 
Reynolds numbers and turbulent -friction coefficients for flat plates 
(ref. 14) are presented in figure 17(a). If transition is assumed to 
occur at the 20-percent longitudinal station of each flight model, the 
friction drag curves would be changed only a negligible amount, since 
the wetted area affected is small compared with the total wetted area. 

The remaining difference in total drag levels appears to be due to the 
effect of Reynolds number on the base drags and experimental accuracy. 

The test Reynolds number (fig. 7 ) of all the models below Mach number 2 
varied between 6 X 10^ and 2 X 10^; and according to reference 15, the 
Reynolds numbers were in a critical range for transition. Reference 15 
also shows that the base pressures and base drags of blunt-base bodies 
may vary significantly at these Reynolds numbers depending on whether 
the boundary layer is fully laminar, fully turbulent, or under transition. 

The base-drag coefficients obtained from the tests and from the semi- 
empirical method of reference 16 for two-dimensional airfoils are compared 
in figure 17(b) • The semi empirical method is restricted to an all- 
turbulent boundary layer at supersonic speeds and requires an analogy 
between the base pressures and peak pressure rise associated with sepa- 
ration of the boundary layer. The comparison indicates that the base 
drag of the model was essentially two dimensional and that the boundary 
layer was turbulent for Mach numbers from about 2 to At Mach 

numbers less than 2 the two-dimensional airfoils have significantly 
higher base drags than the rocket model drag. This result suggests that 
the boundary layer on the model was not fully turbulent at the lower Mach 
numbers and Reynolds numbers. 

Figure 18 presents a comparison of the pressure drags from the rocket 
models, from linearized wing theory, and from supersonic area-rule theory. 
The theoretical pressure drag from wing theory (ref. 17) includes the sum 
of the wing drag with tip effects, fin drag, and interference drag between 
the wing and fins for subsonic leading-edge conditions. The area-rule 
pressure drag was determined in the manner described in reference 18 with 
the assumption that a polyhedron (having the same normal cross section as 
the base) extended downstream from the base to infinity in order to apply 
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the theory. The Fourier series solutions employed for the computations 
of the area-rule drag were convergent. 

The linearized wing-theory pressure drag was in closer agreement 
with the experimental pressure drag than the area-rule drag through most 
of the Mach number range. Although the magnitudes are small, the area- 
rule drag is significantly lower than that from the wing theory. In 
comparison, reference 18 shows nearly perfect agreement for the two 
theories for relatively blunt cones with large bases. 


Coupled Motions 

Typical time-history variations of the aerodynamic coefficients and 
cross plots of C^ against Cy near Mach numbers of 2.3 and 1.8 are 
presented in figures 19 and 20, respectively. The curves defined by the 
data points are generally cyclic but not sinusoidal. It is clearly indi- 
cated that the motions are coupled, the damping in pitch and yaw is very 
poor, and the roll is oscillatory. The cross plots of C^ against Cy 

in figures 19 (g) and 20(g) are not representative of the relative motion 
in pitch and yaw (as they would be for a body having roll symmetry) 
since C N is much higher than C v for the configuration, 
a I |3 

The cross-hatched areas in figure 21 summarize the variations of the 
rolling rates of the model through the Mach number range. The large vari- 
ations in roll-rate amplitude started near M = 2.8 when the first pulse 
rocket disturbed the model. The irregular pattern of positive and nega- 
tive roll oscillations was not changed significantly when the second pulse 
rocket fired. Also shown in figure 21 are the divergence boundaries for 
steady roll, which were determined in the manner described in refer- 
ence 19. It appears that the model rolling velocities are oscillating 
past the unstable regions. 


Simulation of Coupled Motions 

The simulation studies were conducted by use of equations of appen- 
dix B near Mach number 2 by using static stability derivatives from wind- 
tunnel tests and the free-stream conditions of the flight tests. The 
moment derivatives from reference 9 were recomputed about the model 
center of gravity. The pitch damping derivatives and were 

initially calculated from linearized theory (refs. 11 and 12), but these 
values were lowered when they produced too much damping. The damping 
in yaw (^n r an< ^- a -^ so were calculated, in a manner similar to that 

for pitch, with the assumption that all the damping was due to the verti- 
cal fins. The derivatives for were estimated from reference 10, 
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and and were estimated from lifting-line theory. The values 




used for all the stability derivatives are listed in table III. It should 
be noted that two slopes were employed to approximate the nonlinear vari- 
ations of C y and C n with (3, and 


for p, i 
values of 
values of 


C 7 with a. The initial values 
a (3 *Pa 

q were obtained from the flight records; the initial 
p, r, and q were computed from flight data; and, the initial 
a and p were estimated to satisfy the equations of motion. 


L 0 

and 


Consecutive analog runs showed that the simulated motions would not 
repeat exactly for identical initial conditions. The angular velocities 
and their frequencies were similar to those from the test and exhibited 
the type of oscillatory motions shown in figure 22(a) . Small variations 
in the stability derivatives and use of other sets of initial conditions 
did not provide repeatability. The nonrepeatability appears to be due 
to the sensitivity of the equations, with the inputs used, to small com- 
puter errors of the analog. These small computer errors may be compared 
to spurious disturbances during the flight and, in this case, are partly 
responsible for the realistic simulation of the free-f light motions. 

The flight variations of p, r, and q with time are shown in 
figure 22(a) for comparison with three successive analog runs (initial 
conditions in table IV) in figures 22(b), 22(c), and 22(d). The first 
two runs simulate the roll data in magnitude and in the sequence of 
rolling oscillations and reversals. The third run differs mainly over 
the latter half of the run where positive roll oscillations follow the 
second roll reversal. The peak amplitudes and the frequencies of q 
and r were predicted by the equations but their variations with time 
were somewhat different from those of the test. The analog results also 
indicate that the ranges of a and (3 were less than ±0.2 radian 
(about 12°) as the rolling velocities oscillated past the divergence 
boundaries (fig. 17 ) for steady roll rates. 


Several analog runs were made to determine the effect of 

Pa 

on the rolling motions. Oscillatory roll rates were obtained for all 
values between -0.1 and -1.2. Positive values of 0.1 and 0.2 resulted 
in damping of the initial roll rates to very small rates. These runs 
showed also that the rolling motions of the model could not be simulated 
without the C 7 term in the equations of motion. 

1 Pa 

In order to isolate the effect of C 7 on the motions, a simpli- 

fy 

fied three-degree-of-freedom study was made by using equations (b6), 

(B7) , and (B8) of appendix B. All stability derivatives, terms involving 
q and r, and gyroscopic reaction terms were dropped from equations (Bl) 

to (B5) , except C 7 and the angular products |3p and ap. The 

1 3a 
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simplified equations represent a mass, distribution and flight condition 

(near M = 2 in the present case) for which such factors as q S/mV, 

00 

q^Sc/ly, q^Sc^/^VIy, and so on, may be considered negligible with regard 
to the qJBb/ly factor. The angular velocities q and r also were 
assumed to be zero to simplify the equations. 


Figure 25 shows the typical variations of p, a, p, and £ with 
time for the simplified equations using the initial conditions of table V 
and C 1 = -0.82. It can be seen that the rolling oscillations, reversals 


and magnitudes are similar to those of the test in figure 22(a). As in 
the case of the five-degree-of -freedom study, this run also did not repeat 
exactly, showing that the simplified equations with the inputs used were 
sensitive to the small analog errors. Higher initial roll rates resulted 
in repeatable runs with oscillatory roll rates of higher frequency about 
the input value and no roll reversals; lower initial roll rates gave repeat 
able oscillations of lower frequency about zero roll velocity. It appears, 
therefore, that the input conditions which simulated the roll reversals of 
the test were also the critical input values for the simplified equations 
and for the equations used in the five-degree-of -freedom study. Compari- 
sons of a and 3 from the simplified equations (fig. 25 ) with their 
corresponding values from the five-degree-of -freedom study (figs. 22(b) 
to 22(c)) show agreement in their peak amplitudes and frequencies. It 
is evident that the coupled motions, their frequencies and amplitudes, 
were highly dependent on the variation of the effective dihedral deriva- 
tive with angle of attack, the relatively low roll inertia compared with 
the pitching and yawing inertias, and the altitude and velocity of the test 


CONCLUDING REMARKS 


An investigation was conducted to determine the aerodynamic charac- 
teristics of a boost-glide type of hypersonic glider configuration at 
supersonic speeds. The configuration had a 79° clipped delta planform, 
sharp leading edges, a blunt base, and symmetry in two planes. The mass 
distribution was concentrated along the longitudinal axis and the inertia 
in roll was about 4 percent of the pitch or yawing inertias. The test 
covered a free-f light Mach number range from J.k to about 1.2. 

The model was both statically and dynamically stable in pitch and 
yaw at small angles of attack and sideslip near Mach number 5 • When the 
model was pitched to an angle of attack of about 12°, the configuration 
oscillated in pitch, roll, and yaw, indicating coupling of the longitudinal 
and lateral forces. An analog study near Mach number 2 generally simulated 
the magnitudes and frequencies of the angular motions and roll reversals, 
and showed that the rate of change of effective dihedral derivative with 
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angle of attack was a necessary derivative for simulating the coupled 
motions at the altitude and velocity of the test. 

Linearized theory gave reasonable predictions of the static and 
dynamic stability derivatives for small angles of attack and yaw at 
Mach number 3.4. The pressure drag determined from linearized wing 
theory, including interference effects, was in closer agreement with 
the measured pressure drag than the pressure drag from area-rule theory 
through most of the Mach number range. The data indicate that the base 
drag was approximately two dimensional at the higher Mach numbers and 
Reynolds numbers of the test. 

Langley Research Center , 

National Aeronautics and Space Administration, 

Langley Field, Va., September 9, i 960 . 
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APPENDIX A 

ACCELEROMETER CORRECTIONS 


Linear accelerometers which are not mounted on the center of gravity 
measure not only the translatory accelerations but also the accelerations 
due to angular velocities and angular accelerations. Since the linear 
accelerometers in the model were off the center of gravity, it was nec- 
essary to subtract the rotational accelerations from the measured accel- 
erations to obtain the translatory accelerations along the three axes. 

The following equations were used: 

a L,cg = a L,Lo + + r2 ) + y aL , Lo (r ~ Pq) - z a L ^ Lo (4 + P*) 

(Al) 

(A2) 


a Y,cg - ®Y + I 




( r 2 + p2) + Zg^Cp - qr) - 


x a y ^ 


+ pq) 


a N,cg a N,n 


4- i 

g 


Za N,n ^ p2 


q2) _ 


Xa N,n U 


pr) 


ya N ^ n (p 


+ qr) 


(AJ) 


where the distances x, y } and z from the center of gravity to the 
individual accelerometers are given in table II and are identified in 
these equations with subscripts indicating the accelerations measured. 

The corrections for a^ ^ and a^ ^ were made with expressions similar 

to equations (Al) and (A 3 ) . The variations of q, r , and p were meas- 
ured directly from the angular accelerometers and the roll gyro, and 
required no corrections,. The variations of q and r were obtained 
by integration of q and r over several cycles and by assuming that 
they varied symmetrically about zero values. The angular acceleration 
in roll p was obtained by differentiating the variations of p with t. 


The values of a^ and a^ n were used for determining a^ cg 
and a^ cg throughout most of the flight. The values of a^ ^ were 

used only at the higher Mach numbers where the deceleration was greater 
than the range of a^ and a^ ^ served only as a check on c ^. 
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EQUATIONS OF MOTION 


The equations of motion used for the five-degree-of-freedom analog 
study are as follows: 



All gravity terms, as well as certain aerodynamic terms such as C„. 

m (3 

and Cy , have been omitted from equations (Bl) to (B5) . In the calcula- 
x r 4 

tions it was assumed that (3 = -r in equation (B5); thus a slight simpli- 
fication in the last term is possible. All the product-of-inertia terms 
were zero since the principal axis was coincident with the longitudinal 
axis. 


The equations for the simplified analog study with three degrees of 
freedom are as follows: 


4 = -0p (B6) 

0 = a-P (BT) 


P 



(B8) 
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TABLE I 

PHYSICAL CHARACTERISTICS OF ROCKET MODEL 

Wing: 

Planform area, sq ft 5-86 

Base area, sq ft 0.64 

Span, ft 1*77 

Aspect ratio 0. 54- 

Root chord (model length to theoretical tip), ft 5-56 

Mean aerodynamic chord, ft 3*83 

Sweepback angle of leading edge, deg 78-87 

Dihedral, deg 0 

Thickness ratio 0.077 

Taper ratio 0 . 191 

Vertical tail: 

Planform area (2 fins exposed), sq ft O .36 

Base area (2 fins exposed), sq ft 0.05 

Span (total), ft O .85 

Aspect ratio (exposed) O. 5 O 

Sweepback angle of leading edge, deg 76 . 0 

Wedge angle (total), deg 7-67 

Taper ratio (exposed) 0.33 

Mass characteristics: 

Weight, lb 108.0 

Center-of -gravity position, rearward from theoretical tip in 

percent of total length 53-1 

Moment of inertia in pitch, Iy> slug-ft^ 6.93 

Moment of inertia in yaw, I^, slug-ft^ . 7-24 

Moment of inertia in roll, 1^, slug-ft^ 0.27 

Products of inertia, slug-ft^ 0 
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TABLE II 


LOCATION OF INSTRUMENTS 


Instrument 

Distance from 

c.g. 

X, ft 

y, ft 

z, ft 

Longitudinal accelerometer for a L ^ Hi 

0.28 

0.04 

0 

Longitudinal accelerometer for a^^ 0 

.28 

-=f 

0 

1 

0 

Normal accelerometer for a N ^ n 

.075 

0 

0.082 

Normal accelerometer for a^j ^ 

-2.118 

.114 

.059 

Transverse accelerometer for ay 

A31 

1 

b 

Vj4 

0 

Angular accelerometer in pitch for q 

-.144 

-.064 

0 

Angular accelerometer in yaw for r 

.075 

0 

-.071 

Roll-rate gyro for p 

-.406 

! 

0 

0 

Total pressure pickup 

Nose tip 


Base pressure pickup 

Manifolded tube 

on base 


center line 
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TABLE III 

VALUES OF AERODYNAMIC DERIVATIVES USED IN EQUATIONS OF MOTION 


C L ft 

cx 

Cra a • V 

Cm q 

C ra<x 

For 0 < 0 < 0.5 
For 0.05 < (3 . 

C n : 


For 0 < p < 0.1 
For 0.1 < (3 . . 

% 

c i„ > 


For 0 < a < 0.17 
For 0.17 < a . 



1.0 

- 0.12 

-0.05 

- 0.01 


- 0.12 

- 0.18 


0.10 

0.16 

0 


-0.82 

- 0.22 

0.08 

0.01 

-0.06 

0 
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TABLE IV 

INITIAL CONDITIONS FOR EQUATIONS FOR FIVE DEGREES OF FREEDOM 


a 2 .25 

V, ft/sec a 2,100 

Qoot lb/sq ft a 5TO 

a, radians 0.09 

P, radians 0.15 

p, radians/sec 16.75 

q, radians/sec 2.10 

r, radians/sec 0.01 

p, radians/sec^ -282 

q, radians/sec 2 -38.2 

r, radians/sec^ -17-2 


a Average values from 39 seconds to 45 seconds on flight record. 


TABLE V 


INITIAL CONDITIONS FOR SIMPLIFIED EQUATIONS 


a 2.25 

V, ft/sec a 2 , 100 

qoo, lb/ sq ft a 570 

a, radians 0.09 

(3, radians 0.15 

p, radians/sec 12.2 

p, radians/sec 2 -242 


a Average values from 39 seconds to 45 seconds on flight record. 
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Figure 2.- Concluded 
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(a) Planform view of rocket model. 
Figure 3-- Photographs of models. 


L-58-1065a 
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L-58-1064a 

(b) Three-quarter rear view of rocket model. 

Figure 3.- Concluded. 
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Figure 4.- Helium-gun model. 


L-59-489.1 
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L-58-136la 

Figure 5*- Rocket model and boosters on launcher. 
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(a) Free-stream dynamic pressure. 



(b) Altitude. 

Figure 6.- Variations of free-stream dynamic pressure and altitude with 

Mach number for the rocket model. 
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(a) Time interval after model separated from booster. 

Figure 8.- Variations of basic data and free-stream conditions with time 
for three typical time intervals of the flight. 
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(b) Time interval after first pulse rocket fired, 40.8 to 4-1.7 sec 

Figure 8.- Continued. 
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(c) 


Time interval after 


first pulse rocket fired. 


42.2 to 43-1 sec . 


Figure 8.- Concluded. 
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Figure 10.- Variation of C m with C N . M 


3.4. 
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Figure 11.- Amplitudes of q an(3 - Cm as a unction of time. 
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(b) Slope of pitching-moment coefficient. 

Figure 12.- Variations of the static stability derivatives in pitch with 
Mach number at zero angle of attack as determined from wind-tunnel 
tests, flight test, and linearized theory. 
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par radian Cy.i per radian 



(t>) Slope of yawing-moment coefficient. 

Figure 1 6 .- Variations of the static stability derivatives in yaw with 
Mach number at zero angle of sideslip as determined from wind-tunnel 
tests, flight test, and linearized theory. 
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Rocket model 
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(a) Total-drag and computed friction-drag coefficients 


Rocket model iji 

tunnel (ref. 9 ). m 
2- dim. airfoils (ref. 16) Si 


(b) Base-drag coefficients. 

Figure 17.- Variations of total-drag, friction-drag, and base-drag 
coefficients with Mach number at approximately zero angle of 
attack. 
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(a) Total normal-force coefficient. 



(b) Total pitching-moment coefficient. 



40.9 Ui.o l+i.i 4i.2 4 i.3 41.4 


t, sec 

(c) Total axial -force and base axial -force coefficients. 

Figure 19.- Variations of total -force and total -moment coefficients 

with time. M ~ 2.3 . 
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(d) Total side-force coefficient. 



t, B«C 


(e) Total yaving-moment coefficient. 



(f) Total rolling-moment coefficient. 
Figure 19.- Continued. 


CONFIDENTIAL 




• • 


• • •# • 

• • • »• 

• • • • 

• • • * 


»«• 


• • • 

• • • • 

• • • • • 
• • ••• 


CONFIDENTIAL 


• • 
• • • 


• • 


• • 


51 



°Y 


(g) Variation of normal -force coefficient with side-force coefficient. 

Figure 19*- Concluded. 
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(a) Total normal -force coefficient. 




t , sec 


(c) Total axial -force and base axial-force coefficients. 

Figure 20.- Variations of total-force and total-moment coefficients 

with time. M « 1.8. 
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(g) Variation of normal-force coefficient with side-force coefficient. 

Figure 20.- Concluded. 
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Figure 21.- Variations of the amplitudes and ranges of roll rates with 

Mach number for the rocket model. 
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(a) Flight-test motions. Roll rates, angular acceleration in yaw, and 

angular acceleration in pitch. 

Figure 22.- Comparisons of the angular motions obtained from the flight 
tests and from the analog investigation near Mach number 2. 
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